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It is suggested the modification of traditional potential model, in which nontrivial structure 
of inside-hadron vacuum condensate is simulated by geometric properties of inside-hadron space. 
Confinement of quarks is ensured by closed effective (Riemannian) space. Interquark potential 
represents itself Coulomb's low in the effective space (Poisson equation). In the framework of such 
approach the quarks dynamics completely is defined only by metric of the effective space, which in 
turn in conformally-Euclidean case (we consider) is defined by sole phenomenological function. Also 
it is supposed, that inside-hadron vacuum is essential nonperturbative one both on large and on 
small distances and this is taken into account by special parameter in the metric of effective space. 
For final Schrodinger equation the exact analytical solutions for nonrelativistic energy spectrum 
and wave functions of quarkonium are obtained. From the basic principles of geometrized potential 
model and under the perturbation theory the spin-dependent relativistic corrections are calculated. 
Charmonium and bottomonium spectra are simulated. Suggested model gives a very good fit to 
^ ' experimental data: accuracy of spectra fitting makes 4.06 • 10"^ for charmonium and 3.28 ■ fO~^ 

^ I for bottomonium and many predictions are made. On the basis of charmonium and bottomonium 

, analysis conclusions about role of vacuum condensate in hadron structure are done. 

^ : 

^ ■ I. INTRODUCTION. 

> : 

2^ ' Quarkonia are bound states of heavy quarks and antiquarks cc, bb, tt form the special class of objects in the 
, hadrons physics. Experimentally observable quarkonia properties are used for clearing up of static and dynamic 
I ■ quark properties. The quarkonia spectroscopy is studied in detail experimentally j^] and gives in to effective theoretical 
^-H , modeling. The nonrelativistic character of heavy quarks moving within quarkonia allows to apply for calculation of 
^\ • their spectra nonrelativistic Schrodinger equation, in which effective potential of quarks interaction U{r) and quarks 
0^ ' reduced mass as functions of coordinates are selected phenomenologically. Criterion of functions U{r) and /i(r) 
, choice is a good fit of potential model outcomes with experimental data. Smallncss of relativistic effects allows to 
O ^ specialize nonrelativistic potential models by inserting to Hamiltonian effects of spin-orbit and spin-spin interactions 
' , as small corrections computed under perturbations theory. Modern condition of potential models is reflected in the 
review paper ^ . For more detailed information see lectures |^ . 
' Central object in realization of the nonrelativistic potential program is the potential of interquark interaction U{r), 
• • which should take into account effect of absolute quarks confinement within hadron. Its choice is ambiguous in view 
of unsufficient methods of the interquark interactions analysis on large distances. The potential is selected usually 
from some physical reasons in combination with the aesthetic requirements of simplicity. The existing theoretical 



X 



$H indications concerning the form of potential most explicitly are reflected in so-called Cornell potential M 

U{r)^-^+ar, (1) 
r 

where the first term is motivated by perturbative QCD with r — s- with an accuracy to logarithmic corrections, taking 
into account the perturbative vacuum polarization effects; second term is predicted by nonperturbative lattice QCD 
with r —^ oo . Majority of used potentials differ among themselves by character of interpolation on intermediate 
distances, between two limiting conditions reflected in the formula (0). 

Concerning of reduced mass, as function of coordinates, only its asymptotic on small distances, found in perturbative 
QCD, is known. For potential model, however the dependence of reduced masses from coordinates on average and 
large distances is necessary, which unfortunately theoretically is not found. For this reason the reduced mass is 
simulated, as a rule, by constant parameter, which value is determined from the requirement of best agreement of 
theoretical outcomes with experiment. 
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The use of the above-stated procedures allows in whole to achieve a good quantitative agreement of the theory 
with experiment, however, a problem of physical completeness and justification of the potential approach remains 
open. In the framework of the fundamental QCD hadron is an area of reconstructed vacuum - is a cavity inside 
nonperturbative quark-gluon condensate stabilized by valent quarks, residing within it. Moreover, the inside- hadron 
reconstructed vacuum exert influence on quarks dynamics, exactly on reduced mass and interquark potential. Thus, 
quarkonium is complicated self-consistent object, in which quarks condition depends from condition of reconstructed 
nonperturbative vacuum and condition of vacuum inside quarkonium depends from quarks condition. In this situation 
it is worth to suggest the new phenomenological approaches in the theory of quarkonia. The more general potential 
models, at first should contain additional gang of functional and numerical parameters in the interquark potential, 
supposing the interpretation in the terms of nonperturbative vacuum physics; secondly, in these models the problem 
of reduced mass of quarks within nonperturbative inside-hadron vacuum and problem of interquark potential should 
be considered from uniform positions. 

In present work one of the variants of new potential models, possessing by the above-stated properties, is suggested. 
Our approach is based on the consideration of inside-hadron space as manifold with the closed topology. The quarks 
reduced mass and the interquark potential completely are defined by geometric properties of the effective Riemannian 
space; metric of this space represents itself as universal phenomenological functional parameter. We also assume, that 
vacuum inside quarkonium is essential nonperturbative one both on large and on small distances and take this into 
account by specially selected parameter in the metric of effective space, which then occurs in quarks reduced mass 
and interquark potential functional dependencies from coordinates of real (Euclidean) space. 

The paper is constructed as follows. In section II the basic notions about vacuum element in hadron are stated; in 
section III the geometrization of potential model is carried out and its mathematical structure is stated; in section IV 
the generalized Cornell potential, parametric taking into account nonperturbative structure of vacuum on all distances 
is introduced; here the exact analytical solutions for nonrelativistic Schrodinger equation are obtained; in section V 
the relativistic effects are considered; section VI is devoted to the results and conclusions. 



II. THE VACUUM ELEMENT IN THE HADRON MASS. 



In the framework of the qualitative script of the hadron formation the reconstruction of vacuum condensate, 
described quite determined energy, leading to vacuum constituting element in hadron mass, as a necessary condition 
is contained. This vacuum constituting element reads 




{f)dV 



where 

£vacir) = ein{r) - Eout > 0, 

£out is the energy density of outside-hadron vacuum condensate (nonperturbative quark-gluon condensate), Sin{r) is 
the energy density of inside-hadron (reconstructed) vacuum condensate, here the coordinate dependence takes into 
account the spatial heterogeneity of condensate. The mass of hadron is represented as follows: 

M^E,ac. + Eg, (2) 

where Eg is a sum of rest energy, kinetic energy and energy of quarks interactions with each other and with vacuum 
condensate. The birth threshold of hadron agrees (^) is defined not only by quark masses, but also energy Eyac, 
expended to quite defined reconstruction, after which inside a new vacuum the processes of quark excitations births 
with masses become possible. 

III. THE MODEL. 

Let us to start from quarkonium Lagrangian, which in the framework of potential approach reads 

L = l^^{r} (^^^ - C/(f), dll = dr^ + r\d9^ + sin^ 9dip^), (3) 
where fi{f) is the quarks reduced mass and U{f) is the quark-antiquark potential. 
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Having limited by a class of isotropic phenomenological functions /i(r) = m/^(r)/2, U{r), let's rewrite Lagrangian 
(i as 



L=-[jA -Uir), (4) 



df = fir){dr^ + r^ide^ + sin^ edif^)) = -f^pdx'^dx'^ . 



(5) 



The first term in (^) formal coincides with kinetic energy of particle with mass m/2 — const, moving in curved Rie- 
mannian space with metric (H). It means the condensate influence to the quark inert properties can be mathematically 
taken into account by transition from real Euclidean space to effective Riemannian space. Function U{r), appearing 
in (^), within the framework of accepted hypothesis is identified with Coulomb's law in Riemannian space with metric 



^dxf^^^^ dxl^ ~ r^P{r)dr ^ ^ ' dr 



(6) 



From the equation yjm follows 



/x(r) 



m f r dU 
dr 



U{r) = 



dr 



'fir) 



(7) 



The effective potential in (0) is defined with an accuracy to additive constant, which value is parameterized by a limit 
inferior in quadrature. Such model's ambiguity, however, is insignificant, since this additive constant is swallowed by 
vacuum energy ^ , assigning count beginning of energy spectrum. 

For geometrized potential model the hypothesis about additional quark-condensate interaction, which energy pro- 
portional to curvature of the effective Riemannian space is mathematically natural. The insertion this interaction to 
model lead up to Lagrangian 



m dx" dx^ , ^, , 

^ = 4^"^^^^^(^) + '^^('^)' 

where k - phenomenological quark-condensate coupling constant; 

2 „ ,2 



i?(r) 



1 fdf 



/4 \dr 



J3 j,J3 flj. 



(8) 



(9) 



is scalar curvature, evaluated in terms of coordinates (||). 

For Lagrangian the corresponding Schrodinger equation is 

-— 7"''VaV;3 + Uir) - KRir)] * 
m / 

where Vq is covariant derivation operator in space with metric (^. After variable separation 
equation (hoh is reduced to equation for the radial part of wave function 



(10) 



(11) 



1 1 d 2r dRnL 

-r f- 



m r"^ p dr '' dr 
with invariant normalization 



1 L{L + l) 
m r'^p 



dr 



2k 



1 



J3 J.J3 flj. 



RnL — EnhRnL (12) 



^/ld' 



RlL{r)f{rydr = l. 



(13) 
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From the view point of physical representations, explained above, the choice of non-Euclidean measure of integra- 
tion in (^3|) is finishing mathematical operation under the account of condensate influence on spatial localization of 
quarkonium state. 

The radial coordinate r, used in (||) - ([T^), is common both for the real Euclidean space and for the effective 
Riemannian space. Therefore, arguing of problem about radial dependence of quark-antiquark potential is carried 
out in terms of this coordinate. However, in mathematical research of the equation (|l^) it is more convenient to use 
dimensionless coordinate p and phenomenological function 



f{r)dr , x{p) = -r{p)f [r{p)] , Pmax = p(oo) 
a 



(14) 



The metric (|^), curvature and Coulomb's law (Q) become: 



dl'^ = {dp^ + x^{p) {de^ + sin- 



Rip) = ^ 



)] = japdx'^dx'^ 
p 

, U{p)^9^ 



dp 



x^{p) 



(15) 



In ( p^ and in the further primes derivative are designated on p. The physical and the geometrical meaning of 
dimensional parameter a become clear after giving representations for potential, satisfying to condition of quarks 
confinement within hadron. 

Further, it is convenient to turn to new radial function with simple normalization: 



R„ 



1 



a3/2x(p) 



4>nL{p) 



(l>nL{p)dp = 1 



(16) 



and to dimensionless quantities: 



-En 



27 



g^ma 



2 Km 



(17) 



Now equation (|T^) corresponds as follows: 

p 



EnL - 27 



dp L{L + 1) 



X 



X 



X 



Po 



/2 



(18) 



IV. THE GENERALIZED CORNELL POTENTIAL AND EXACT SOLUTIONS OF THE SCHRODINGER 

EQUATION. 



In quarkonium geometrized potential model the radial dependence of the effective potential is set simultaneously 
with the metric and the curvature of effective Riemannian space by sole phenomenological function: 



fir) 



r2 dU 
g'^ dr 



(19) 



Among of potentials satisfying to condition of quarks confinement within hadron the potentials of whirlpool form 
stands out of the simplicity: 



U{t) 



2a 



l/k 



1/k 



(20) 



Where 2a is characteristic quarkonium size ("diameter"); k is positively defined numerical parameter. With k ~ 1 
expression (EOh passes to Cornell potential [0, joining known asymptotics. Euclidean-Coulomb asymptotic of Cornell 
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potential with r — !■ corresponds to supposition that properties of inside-hadron vacuum in a neighborhood of point 
r = hmiting are closed to properties of perturbative vacuum. It is possible to expect, with k ^ 1 generalized 
Cornell potential (|2^) simulates situation with nonperturbative inside-hadron vacuum condensate. For potential ( |20| ) 
the evaluations under the formulas (^9|), (^4|), (|l^) give: 



f{r) 



k 



' 2a \ k' 



' J_\~ 

^2a) 



r = 2atan''-^, < p < irk'^ 



X{p)^ksm—, U{p) 



sin^^ 



(21) 



The equation (|18|) and normalization (16) become: 



+ I k^er^L + 3A + 1 + 2"ik'^ cot 7? \y- \k^L{L + I) + + (j)„L = 0, 

t sm ?7 J 



(22) 







1, 



(23) 



where 

Tj = p/k'^, lyiiy + l) = X{1- k^). 
The solutions of the equation ( p^ ) satisfy to condition of absolute quarks confinement within hadron, if: 

1 /I 



X{1 - k')iy = -- + J - + X{1 - k^) >0 



(24) 



The inequality ( p4| ) is the sole restriction on the phenomenological parameters of geometrized model appropriating to 
generalized Cornell potential (po|). With k — 1 the effective geometry, as it is visible from (|2l| ), represents a closed 
homogeneous and isotropic space with constant positive curvature R — The equation (|22|) in this case describes 

hydrogen-like system in such space: 



a" 

Kl 



£nL + 3A + 1 + 27 cot p - 



L{L+1) 
sin^ p 



= 0. 



(25) 



From symmetry reasons is beforehand obvious, that spectrum of such system depends only from the sum of radial 
and orbital numbers N = n + L, i.e. degenerates on orbital number. 

With k ^ 1 the degeneracy is removed, but mathematical structure of the equation for wave function on comparison 
with (^Bf) practically does not vary. Really, after insertion variable 77 and intermediate auxiliary parameter I: 



l{l + 1) = v{v + 1) + k^L{L + 1) , 1 = — 




the equation (E2[) corresponds as follows 



k'^£nL + 3A + 1 + a-fk^ cot 77 



^(^ + 1) 
sin^ 77 



k^L{L + l), 



= 0, 



(26) 



and differs from (E3) only by overdetermination of dimensionless quantities 7 ^k^, SnL £nLk'^ , L I. The 



mathematical uniformity of the equations (|2^) and (^6|) does not cancel, however, their physical distinctions: in ( |26| ) 
I is not orbital quantum number, accepting only integer value. In model with potential (|20| ) and with additional 
quark-condensate interaction, proportional to curvature, with k ^ 1, parameter I with realization of the condition 
(24) can accept any positive values even for states with real orbital moment L — 0. The solutions of the equation 
(26) are energy spectrum: 
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and exact eigenfunctions 



n 



d {I + s) cos 77 
di] sin 77 Z + s 



(27) 



(28) 



In (|2^) line-ups of operators is ordered in such a way, that operators, corresponding to smaUer values s, stand to the 
left of operators, corresponding to large s. The auxiliary functions tpnL satisfy to equations 



d {l + n)cosr] 
di] sin r] I + n 



VnL = 0, 



having exact solutions. 



IfinL = CnL exp(- 



I + n 



) sin"+' 77 



where CnL - normalizing constants, computed numerically from normalization (p3|). 

Turning to dimensional quantities under the formulas ( p7| ) and choosing additive constant in the correspondence 
with (0), we shall receive from (uTh the expression for nonrelativistic quarkonium spectrum: 



Mni 



rriQ + mi n — 



fc2L(L + 1) 



1712 



{n~\ + ^{. + \f + k^L{L + l)^' 



(29) 



where parameters are defined as follows: 

mo — 2m + 27r'^a'^fc^£o, 'Tii — 



2 3; 4 1 "^47 



4;,4 



Energy, expended to reconstruction of nonperturbative vacuum condensate inside quarkonium, here is presented as 
Eyac — 27r^a'^fc^eoi where V = 27r^a^k^ is quarkonium volume, calculated under the metric of effective Riemannian 
space; £0 - parameter quantitatively describing vacuum reconstruction. 



V. FINE AND HYPERFINE SPLITTING. 



For the analysis of the relativistic effects it is necessary beforehand to agree about interpretation of nonrelativistic 
levels ( p9| ) . There are two possible interpretation. According to the first of them the eigenvalues of the nonrelativistic 
Hamiltonian are centers of gravity of the spin triplets i.e. ortho(5 = l)quarkonium levels, according to the second 
the nonrelativistic Hamiltonian gives the centers of gravity of the full spin multiplets. For the states with L ^ 0, this 
interpretations are equivalent and the difference occurs only for the S'-states. In the framework of potential ideology 
there is not justification for one interpretation rather than the other. From the practical reasons that the masses 
of the spin singlets for bottomonium are not known we interpret nonrelativistic levels as the centers of gravity of 
ortho quarkonium . 

In present work we shall be limited by consideration only spin-dependent interactions - Hsd- Following to standard 
procedure, based on Breit-Fermi Hamiltonian Q and in conformity with general principles of geometrized model, 
explained in Section HI, let's write Hsb in the covariant form: 

Hsd = -^e„^7(^Q + (V'^C/) (V) - 1-^{s'^qV ^){s^-V p)U + ^ [V"V„C/] (7./3Sq4) , (30) 

where Vq is the covariant derivation operator; Sq = ^ctq is (anti)quark spin operator. The matrixes ctq are defined 
in the effective Riemannian space by condition 

^q4 + '^q^q = '^i"^ (31) 
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and covariant derivative^ 



V/3a§ = 

Representing the metric (|[) as follows 

df = f{x, y, z) [dx^ + dy"^ + dz^] = -fa.pdx^'dxP 



(32) 



and further taking into account the spherical symmetry of potential and ratio (19) Haniiltonian (^ 

HsD = Hls + Ht + Hhyp, 

where 



is reduced to 

(33) 

(34) 



Ft 



1 



2m? \r^f 



\r_df_ 
3fdr 



iS'^x^XS^xp) 



~ S Sq 



(35) 



6{r) 



(36) 



here 5" = Sq + s^. Hls ( P^ ) and i?T ( |35| ) are accordingly the spin-orbit and the tensor interactions, responsible 

for the fine structure of orthoquarkonium levels, Hhyp (|3^ ) is the hyperfine interaction, splitting the states ortho and 
para(5' = O)quarkonium. The number ( |36| ) corresponds to our assumption that the nonrelativistic Hamiltonian 

gives the orthoquarkonium levels. 

The numerical values for the fine and hyperfine splitting of nonrelativistic levels are found by calculation of matrix 
elements for the operators Hls, Ht, Hhyp on wave functions of corresponding states and by virtue of factorization 
( |ll| ) the determination of the matrix elements of spin and spatial operators can be carried out separately. 

Let us to consider at first the fine splitting. The matrix elements for the spin operators in (p4) and (BSf) are 



K 



(LS) _ 
JL — 



(J,L,1|LS|J,L,1) = -[j2-L2-2] 



(37) 



J,L,l ) = -- 



6 (LS)-^ + 3 (LS) - 4L2 

4L2 -3 



(38) 



where = Xa/r, J = L + S. After transition in (34) and ( |35| ) to coordinate rj, we can write the final expressions 
for numerical corrections as follows 



XfpiLS) _ . ly-iLS) A _ 



1 



8^ J sinr/tan2^-f '^"^ 



dij 



(39) 



^-^iiJ = BulK^^l, BnL 



(T) 



3ff2fc 



sin 77 tan^*^ ^ 



2fc + cos?7\ ,2 , 



(40) 



Let us to note that expressions (|39|) and ( [401) do not contain any free parameters. 

When analyzing hyperfine interaction, we suppose additional quark-condensate correlation, conditioned by influence 
of spin-spin states of quarks on the state of nonperturbative vacuum inside quarkonium. We understand that is 



^Expressions (pll) and (p2|) are nonrelativistic limit of fundamental spatial-spin connection 7*7'' -I- 7'"7' ~ 2g''' and Vijk = 
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reflected in the form of interquark potential, formed with active participation of this vacuum. For account with 
supposition we replace phenomenological parameter k (appropriating to orthoquarkonium) in the potential (^) by 
a new phenomenological parameter ko appropriating to paraquarkonium, which value is fixed from the adoption of 
the experimental data. Since there is not experimental data for parabottomonium we shall not consider the spin-spin 
interactions in this quarkonium. 

Further we shall notes, that the S function in ( |36| ) comes from the unnatural nonrelativistic reduction and will 
become a smooth function with a finite range if this is calculated correctly. The standard way to calculate Hhyp is to 
replace 5{f) by a smeared function ||^, which we assume 



Sir) 



sin {fir} 



lim S{r) = d{r) 



where fi is the range of smearing, which is a free phenomenological parameter, we choose M = ^5^- 
In terms of coordinate rj the numerical values of hyperfine corrections are 



SE. 



(hyp) 



2g^ /•sin{itan'=(i)} 



irn^a^k 



SVC? T] 



(41) 



The criterion of adaptability of such mode calculation Hhyp can be following. For a P state, the center of gravity 

(COG) of triplet l'^P(cOG) must coincide with singlet l^Pi if Hhyp is pointlike interaction, since |4'(0)r = 0. From 
experimental data for charmonium, one finds 



l^P 



(COG) 



Our calculation under the formula (^Tj) give 



l^P 



(COG) 



l^Pi = 1 Mev 



l^Pi = 5 Mev 



This is a reasonable accuracy and this means the range of Hhyp is reasonable. 

Finally, the paraquarkonium levels are determined in two stage: on the fist stage the nonrelativistic levels 
parameter fcg are found and on the second the corrections (pl| ) are calculated. 

Wave functions (/)„i used in expressions (|39|), (EO) and (|41|) in correspondence with (pi) have a following form 



) with 



kh = CiL exp 



jk*r] 

'ITT 



sin'"*"^ ri; 



4'2L = C2L- 



exp 



(-^) 



• 1+2 
sm 77 



(/ + 2){l + l)smr] 



[{21 + 3)7A:^ sin 77 - (2^^ + 91^ + 131 + 6) cos ?]] ; 



t>3L = C3L 



exp 



(-^) 



sin'+^ 77 



(Z + 3)2(Z + l)sin^?7 



{Yk'^l + 107^A:'*45 - 931 - 651^ - 2r - 191^) sin^ r? 



+ (4?^ + 44/4 _^ j^gy^a _ 3gg^ _^ 3g]^p _^ ^^35) ^^^2 ^ _ (8^fc4;3 _^ ^Q^^^f _^ 1267^'* + 907/0^) cos 77 sin 77 
and are normalized according to (E3) . 



VI. RESULTS AND CONCLUSION. 



In Table 1 the values of parameters, appearing in expression for nonrelativistic spectra (29) are presented, which 
were used at all stages of charmonium and bottomonium spectra numerical modeling. In Table 2 the values of 
quarkonia physical parameters are presented. Table 3 contain outcomes of charmonium and bottomonium spectra 
modeling for experimentally observable levels. On Fig. 1, 2 in standard representation the complete charmonium and 
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bottomonium spectra, including experimentally known levels and our theoretical prediction for not yet opened levels 
are presented. The levels have spectroscopic designation ^^'^^Lj and are located in corresponding to experimentally 
fixed quantum numbers J^*-^. 

Accuracy of spectra theoretical modeling was estimated under the formula 



\ 



N 

E 



M, 



(th) 



M, 



i{ex) 



M, 



i(ex) 



Ml. 



(42) 



where N - number of quarkonium levels, on which is conducted comparison theories Mj(j^) and experiment Afj(e^) 
behind elimination of ground state of orthocharmonium Afis = 3097 MeV and orthobottomonium Mis — 9460 MeV, 
accepted for a count beginning. For all experimentally known levels of charmonium TV = 11; in this case the calculation 
under the formula (42) give 6['j^ = 4.06 • 10~^. Levels below charm threshold is estimated under the same formula 
(El) with N — 7; here we have s!f ^ — 5.05 • 10^^. For all experimentally known levels of bottomonium N — 11, here 

= 3.28 • 10"^ For levels below bottom threshold = 8, ^''^ = 0.96 • 10"^ In Table 4 the accuracy of spectra 
modeling of the several best models known for us, calculated under the formula (|4^), are presented. 

From numerical values of charmonium and bottomonium parameters, presented in Tables 1 and 2, follows the 
quarkonia are essentially nonperturbative objects. Foundations for this conclusion serve, at first large value of a 
coupling constant on quarkonium scale; secondly, essential difference of parameter k from one, showing the absence 
of Coulomb's asymptotics with r — > 0, that, in turn, means essential nonperturbative properties of vacuum inside 
quarkonium on all scales. The last conclusion is confirmed by comparison of parameter Eq, describing a power of 
vacuum condensate partial reconstruction inside quarkonium, with the module of outside-vacuum condensate energy 
density §. 



0.3- (0 



Ga /^anu 



0) w {27bMeVy 



As we see, the deconfinement of quarks within quarkonium take place already with eo/|£out| ~ 0.03 for charmonium 
and with eo/|eo«t| ~ 0.05 for bottomonium. Conclusion about more force reconstruction of nonperturbative vacuum 
inside bottomonium in comparison with charmonium is qualitatively coordinated with values of nonperturbative 
coupling constant g^: in bottomonium it is less, than in charmonium. 

In the end let us to make conclusion: in the framework of our model we can quantitatively see, that nonperturbative 
vacuum play key role in hadron structure, it is necessary and defining ingredient of quarks dynamics and, therefore 
observable properties of quarkonia. By justification of this conclusion the good fit to experimental data serve. 
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TABLE I. Parameters of model. 



Parameter 


Charmonium 


Bottomonium 


V 


0.158 


0.139 


2fco 


1.291 




21? 


1.179 


1.360 


mo 


3628.209 MeV 


9963.116 MeV 


nil 


49.348 Mev 


48.136 Mev 


m2 


800.825 Mev 


734.050 Mev 



TABLE II. Parameters of quarkonia. 



Parameter 


Charmonium 


Bottomonium 


rn 


1700 MeV 


4905 MeV 




2.328 


1.138 


a 


5.856 ■ 10"^ MeV"^ 


3.028 ■ 10"^ MeV"^ 


eo 


(113.442 MeV)'' 


(129.298 MeV)" 


A 


57.243 MeV 


68.278 MeV 



TABLE III. Fitted levels of quarkonia (MeV). 



Charmonium 



Bottomonium 



State 



Data 1 1 1 



This work 



State 



Data 



This work 



l^Si 

I^PI 

l^Po 
l^Pi 



2'So 
2^81 
l^Di 
3=^81 
2=^Di 
4=^81 



2980 

3097 

3526 

3415 

3511 

3556 

3594 

3686 

3770 

4040 
4159 ± 20 
4415 ± 6 



2974 
3097 
3512 
3413 
3517 
3556 
3647 
3686 
3772 
4040 
4127 
4435 



r8i 

l^Po 
l^Pi 

1='P2 

2^81 
2^Po 
2^Pi 
23P2 
3'8i 
2^Di 
4=^81 
3^Di 



9460 
9860 
9892 
9913 
10023 
10232 
10255 
10269 
10355 
10580 
10865 
11020 



9460 
9860 
9901 
9916 
10023 
10222 
10254 
10266 
10363 
10522 
10745 
10936 



TABLE IV. Accuracy of spectra modeling (%). 



Accuracy 


EQ i 


GJ H 


MZ 


ZOR 


BBZ 


This work 


"7 


6.36 


1.9 


2.37 


9.83 


5.34 


5.05 


"11 








9.83 




4.06 




5.15 




0.14 


3.25 


2.04 


0.96 
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